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OPTIMAL N-BURN MULTIORBIT INJECTION 


by 

Robert L. Singleton, Jr. 
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degree of Master of Science. 

ABSTRACT 

A method of minimizing gravity losses in orbital escape 
maneuvers, known as "perigee propulsion" or "multiorbit injection", 
has been found to have a useful application in transfer problems 
involving high specific- impulse powerplants such as nuclear rockets. 
Unlike conventional chemical rockets which typically have thrust 
phases that can be approximated by velocity impulses, new- genera- 
tion powerplants require powered phases of finite duration to 
effect a required velocity change, and the impulsive approximation 
is no longer valid for evaluating the performance of the vehicle. 

A theory has recently been proposed, however, which establishes , 
ground rules for computing characteristic velocity losses over an 
optimally- steered, finite- thrust trajectory on the basis of the 
velocity impulse required for the same maneuver. Utilizing this 
theory, this thesis develops a procedure for computing gravity 
losses over an N-burn multiorbit escape trajectory of specified 
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final energy and presents a technique which can be readily and ef- 
ficiently employed to predict optimal burn schedules for time- open 
and time-fixed multiorbit escape maneuvers. 
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CHAPTER I 


INTRODUCTION 


An interplanetary orbital transfer, such as injection from 
an orbit about the earth into a hyperbolic trajectory to Mars, 
can theoretically be effected without gravity losses by an impulsive 
velocity change at perigee. Using conventional chemical- fuel 
propulsion, the velocity change required by the transfer, i. e. , 
the difference in periapsis velocity at the junction of the two orbits, 
can be achieved with a high- thrust burn, the duration of which 
is typically much less than the period of the initial orbit. Hence, 
for many orbital transfer problems, a single velocity impulse is 
very often a good approximation for the added velocity increment. 
The impulsive approximation is not valid, however, for vehicles 
requiring finite burns to effect a required velocity change, where, 
by "finite", it is presumed that the burn time is not negligible 
with respect to the period of the orbit and that the thrust segment 
is distributed over a significant portion of the flight path. Finite 
thrust characteristics must be taken into consideration in missions 
involving new- generation powerplants, such as solid- core nuclear 
rockets or electric propulsion devices. For example, a typical 
low- acceleration nuclear- powered injection maneuver using 
continuous thrusting would require a long, spiralling trajectory 
with much of the energy being applied in regions of low velocity. 

A critical factor in the determination of the initial thrust- 
to- weight ratio is gravity losses, which are defined to be the 
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difference in the impulsive velocity change AV^ required for an 
injection maneuver of specified excess hyperbolic velocity and 
the characteristic velocity AV-p required by the vehicle to perform 
a maneuver of the same magnitude using finite thrust. Gravity 
losses, therefore, are interchangeably referred to as character- 
istic velocity losses AV*, given by 

AV* = AV f -AV : (1-1) 

Energy addition efficiency is defined to be maximum for a given 
transfer when there are no characteristic velocity losses and 
when energy is applied in the region of highest velocity, as in 
the case of an impulsive transfer at perigee. Since the losses 
accumulated during high- thrust burns of short duration are 
very small, energy addition efficiency is near- optimal, and 
one- impulse transfers using conventional powerplants are a 
viable means of interplanetary maneuvering. For new generation 
powerplants requiring finite thrust, however, an alternate thrust 
program to minimize gravity losses must be proposed if such 
devices are to be applied economically and efficiently to inter- 
planetary missions. 

A thrust program which dramatically reduces gravity 
losses over finite- thrust transfers has been proposed under names 
such as "perigee propulsion" [9] and "multiorbit injectiori'[ 5] . 
These multiburn- multiorbit techniques maximize energy addition 
efficiency by applying thrust intermittently in regions of high 
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velocity (e. g. , at perigee) and allowing the vehicle to coast under 
the influence of gravity between burns. The resulting flight 
path, shown in Figure 1, consists of a series of coaxial ellipses 
of sequentially increasing energy, along each of which the vehicle 
falls until it reaches the desired position with respect to the 
next perigee, where thrusting resumes and transfer to a higher 
orbit is effected. " When the vehicle arrives at perigee with 
sufficient energy to effect the required injection with a short burn, 
the multiorbit injection maneuver is complete. 

The advantage of using the multiorbit technique for the 
injection maneuver is clear. Whereas a continuous thrust, 
spiralling trajectory typically requires initial thrust- to- weight 
ratios near 0. 5, the multiorbit injection maneuver can make 
use of accelerations less than 0. 1 g [9]. Furthermore, a speci- 
fied payload can be powered by a smaller engine using a multiburn 
thrust schedule, or the payload for a given mission can be increased. 
This improvement in energy addition efficiency is gained, however, 
only at the expense of longer required transfer times. Neverthe- 
less, for interplanetary missions requiring new generation power- 
plants, transfer times on the order of several days are still very 
small when compared to overall mission times measured in 
months or years. 

The purpose of this study is to develop a simple method 

for arriving at the optimal multiorbit burn schedule for both 

"Pericenter is not at the same position for each orbit as is 
implied in Figure 1. Each burn is accompanied by a small radial 
displacement as described in Chapter 2. 
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Figure 1.1 Five-Burn Multiorbit Injection 


(10^ mi ) 
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time- open and time-fixed injection maneuvers from a circular near- 
Earth parking orbit to a hyperbolic escape trajectory of specified 
excess hyperbolic velocity V . The optimal burn schedule is 
defined to be that thrust program which minimizes the gravity 
losses for a multiorbit injection maneuver of a specified number 
of burns, subject to the initial orbital parameters, the final 
time and energy constraints, and the nature of the characteristic 
velocity approximation. In the following analysis, the effects 
of a third body will be ignored, and a simple planar model for an 
orbital transfer in an inverse- square central field will be utilized. 

The following chapter will develop the equations used to 
calculate gravity losses over an N-burn multiorbit trajectory and 
will outline the assumptions required to permit a closed- form 
computational technique. Also included in Chapter 2 are the im- 
pulsive velocity and primer vector solutions for the injection 
maneuver. Chapter 3 contains a derivation of the approximate 
solution to the optimal time- open transfer and a presentation of 
an efficient algorithm for generating optimal N-burn trajectories 
for both time- open and time-fixed injection maneuvers. Chapter 4 
presents a description of the procedure used and a summary of 
the results obtained for a typical multiburn injection maneuver. 

A set of conclusions are presented in Chapter 5. The Appendix 
contains a derivation of Robbins' "impulsive approximation", 
a listing of the orbital and powerplant data for a typical injection 
maneuver, and an outline of the computer simulations used in 
this study. 15 



CHAPTER II 


GRAVITY LOSSES ON A MULTIORBIT 
INJECTION MANEUVER 


2. 1 The "impulsive Approximation" for the Characteristic Velocity- 

Earlier studies [5, 6, 9] have considered the multiorbit 
injection problem and have computed characteristic velocity losses 
using non-optimal steering criteria and numerical integration 
over the thrust trajectory. This analysis, however, will develop 
a procedure to obtain the optimal N-burn thrust program which 
incorporates not only optimal steering but. also a convenient and 
simple closed- form method for calculating gravity losses. The 
groundwork for this analysis has been laid out in a report by 
H. Robbins [12], in which a method is presented for predicting 
characteristic velocity losses over a finite- thrust trajectory on 
the basis of the velocity impulse required to effect the equivalent 
maneuver. Robbins' approximation is derived in Appendix A. 

The result of Robbins' approximation, as derived in the 

- 1 < 

Appendix, expresses the characteristic velocity loss AV" over 
a segmented finite- thrust trajectory in terms of the primer 
vector solution and the velocity impulses over the equivalent 
impulsive trajectory. The steering is assumed to vary linearly 
over each burn, and the centroid time for each burn corresponds 
to the time of occurrence of the equivalent velocity impulse. 

The result is given by 
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n 

AV* =1 I (- X * G X - |ij 2 ) M 2 (2.1-1) 

2 k=l ~ “ ^k 


where \ is the primer vector of Lawden [10], is the second 

— k 

moment of thrust about its centroid, and n is the number of burns. 

The second moment can be calculated from the expression 



_ 1 _ 

12 


AV i T b 2 
D k 


(2. 1-2) 


where AV t is the equivalent velocity impulse on the k th burn 
X k 

and T, is the duration of the k th burn. (2. 1-2) is exact for 
D k 

constant thrust over the burn segment; however, to account for 
non- constant acceleration during the burn, the above equation 
can be multiplied by the correction factor f, which for the k th 
burn is given by 


f 


k 


6c 

AVy 

x k 



+ coth 




= 1 


60 




(2. 1-3) 


where c is the exhaust velocity. 

Although the characteristic velocity on the k th burn 

is calculated as a function of the equivalent velocity impulse 

AV for that burn, it is important to note that the k th finite- 
X k 

thrust segment is actually centered about a generalized velocity 
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impulse AV^ , 
G k 

impulse times 


which has position discontinuities at the 
given by- 


d k = -\(t k )M 2 (2.1-4) 

k 

and which differs in magnitude from AV T as follows: 

Tt 

AV q - AV = X(t k ) ' d k = - |X(t k )[ 2 M 2 (2.1-5) 

k k k 

Since X (t k ) is perpendicular to X (t k ), the position discontinuity 
is radial at each perigee burn and its magnitude is of the order 
of |X (t k ) [ . Thus, the displacement of one perigee with respect 
to another on an optimally steered and timed multiorbit trajectory 
is small enough with respect to the semimajor axis of the orbit 
to be neglected. 

It can be shown that the first term in parentheses on the 
right hand side of (2. 1- 1) is equivalent to the gravity gradient 
-g , where x is the mean direction of thrust acceleration, that 

AA 

* 2 

is, along_X Ct k ). The second term, jx| , is equal to the square 
of the optimal turning rate of the thrust vector. For short burns, 
a constant average value for g xx is a good approximation; and, 
for the simple orbital model described, this term is given by 


-g- 


xx 


= 4 - (1 “ 3 sin2 e > 

r d 


(2. 1-6) 
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where jx is the gravitational constant, r is the distance from 
the center of attraction, and 9 is the angle between the mean 
thrust direction and the local horizontal. Therefore, for a 
n-burn multiorbit trajectory in a planar, inverse- square central 
field, Robbins' approximation yields the following equation 
for computing the characteristic velocity loss: 


n 


AV = 


24 


1 [Jjpd- SsinV |x| 2 ]aV f k T b 
1 r & 


k= 1 


(2.1- 7) 


2. 2 Thrust and Orbital Constraints 


The multiorbit injection maneuver considered in this study 

consists of an N-burn transfer from a circular orbit of radius 

R (close to the earth) to a hyperbolic escape trajectory with 

specified V^. The transfer profile describes a series of (N-l) 

intermediate elliptical coasting orbits between the initial and 

final burns. The transfer time T is equal to the sum of the 

s 

periods of the (N-l) coasting orbits. Each perigee burn is 
executed such that its centroid occurs at the junction of the 
circular orbit and the (approximately equal) perigees of the 
N succeeding conic orbits, and the thrust acceleration vector 
is applied parallel to the primer vector, such that optimal 
timing and steering requirements, respectively,, are satisfied. 
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Having specified N, R, and , sufficient information is 
available to compute the characteristic velocity loss for a multi- 
orbit escape maneuver from an Earth orbit. Certain constraints 
must be imposed on the analysis, however, to assure validity 
of Bobbins' approximation and a closed-form calculation of the 
performance penalty. First, r and 0 are assumed to remain 
constant over each burn. This assumption places no further 
restrictions on the problem, since it is consistent with Robbins' 
short burn criterion. Therefore, assuming that the position 
discontinuities d^ are much smaller than the radius of the initial 
orbit, and specifying r = R and 0 = 0 for each perigee burn, the 
characteristic velocity loss for an N-burn multiorbit injection 
maneuver is given by 





f k T b, 


( 2 . 2 - 1 ) 


Exactly what is meant by a "short" burn in the analytical sense 
is not given in [12]. However, Bobbins' analysis demonstrates 
that, in order for a burn to be short, the dimensionless quantity 
co g must be sufficiently small, where to , the Schuler Frequency 
is defined as 


to 


s 



( 2 . 2 - 2 ) 


Furthermore, by a rough analysis. Bobbins finds that the short 
burn approximation is quite good for values of to T, as large as 

’.iv S D, 
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unity and progressively worse for values beyond this arbitrary 
value. One indication of the range of flexibility offered by this 
semi- arbitrary short-burn constraint is the maximum allowable 
burn time permitted by it. Comparing the maximum burn time 
given by 


T 

b max 




with the period of the circular orbit. 


(2.2-3) 


T o = 27T ^Z~ 

The following limit is obtained for a short burn:- 

T b< - 16T o <2 - 2 - 4) 

k 

In other words, using the short burn criterion based on to , 

s 

Robbins' approximation is valid on thrust trajectories over arcs 
up to a magnitude equivalent in time to a trajectory covering 16% 
of the circular orbit. This result is rather optimistic, if not 
surprising. Nonetheless, Robbins' short burn criterion has been 
shown to give accurate results and will be used as an upper limit 
in the optimal N-burn multiorbit injection analysis. 


For a time-open transfer, that is, one for which the total time 


between the initial and final burns is unspecified, the orbital 
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constraints are the requirements that (1) the velocity after the 
last burn will satisfy the final energy constraint and (2) the 
velocity on the next-to-last burn (k = N - 1) will be less than the 
escape velocity, given by 

V E = V^Vq (2.2-5) 

where V q is the velocity of the vehicle in circular parking orbit 
given by 


V 

o 



( 2 . 2 - 6 ) 


The first requirement will be described in the next section. The 
latter criterion is necessary to constrain to transfer to the 
specified number of burns N. 


2. 3 The Velocity Impulse and Primer Vector Solutions 

For an impulsive transfer from a circular orbit to a hyperbolic 
trajectory, the required velocity impulse AV^ is equal to the 
difference in the periapsis velocity of the hyperbolic trajectory 

and the initial circular velocity V . The energy equation 
for a hyperbolic orbit is given by 



(2.3-1) 
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where a, the semimajor axis is taken to be negative. If the 
escape velocity at infinity (r — > <») is specified to be V , a is 
determined uniquely, and the hyperbolic velocity at perigee 
(r = B) is given by 


or 



7T 


2 

+ V 

CO 


(2.3-2) 





(2.3-3) 


The total required impulsive velocity can then be written as 


AV 


I 



V 

o 



+ V 



(2.3-4) 


The primer vector solution for this transfer can be obtained 
in a simple manner from an examination of the necessary 
conditions along an optimal thrust trajectory. In order for energy 
addition efficiency to be 'maximized along an impulsive trajectory, ' 
the velocity impulse must occur at a local maximum of the primer 
vector time history, and, immediately after the impulse, necessary 
conditions for the trajectory to be optimal are given by the 
relations [10] 

X a v 

X«g (2.3-5) 
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where, immediately after the impulse, the magnitude of the 

2 

gravity vector is (ju/R ) and the velocity is the periapsis velocity 
of the hyperbola Vjj . From the conditions given by (2,3-5), 
the following relationship can be written: 

k £ 

X V (2.3-6) 


Thus, since the magnitude of the primer vector immediately after 
the impulse is defined to be unity along an optimal trajectory, 
(2.3-6) can be written as 



(2.3-7) 


and the primer vector solution to be used in the calculation of the 
performance penalty (2.2-1) is given by 

(2.3-8) 
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CHAPTER III 


OPTIMAL N-BURN THRUST PROGRAMS 


3.1 Approximate Time -Open Solution 


The optimal burn schedule for the time-open multiorbit 
injection maneuver is that thrust program which, subject to the 
appropriate constraints, minimizes the characteristic velocity- 
loss over the entire trajectory. Using the primer vector solution 
obtained in Section 2. 3, the characteristic velocity loss given by 
(2. 2-1) can be rewritten as 


AV 




(3.1-1) 


which can be reduced to 


AV 


N 

C -I 

k=l 


f T 
I, k b, 
k k 


(3.1-2) 


where the constant C ro , specified by the transfer constraints, is 
given by 


v «, E + M 


24 R^ \v 2 R +2 n 


(3.1-3) 


For a single burn, the burn time T^ for a constant thrust 
powerplant can be expressed in terms of the equivalent velocity 




change AV over the burn segment, i. e. , 


T, 



( 3 . 1 - 4 ) 


where m is the initial mass of the vehicle, c is the exhaust 
o 

velocity, and F is the thrust magnitude. For a multiorbit- 

multiburn trajectory, the burn time for the kth powered phase 

can be expressed by taking the difference in the burn time befope 

and after the burn. If T. and T, are defined to be the times 

b i b i+i 

at the beginning and end, respectively, of the kth powered phase 

(i = k), and if AV. and AV\ + ^ are the corresponding values of the 

total impulse up to times T, and T, , respectively, then the 

D i l+l 

velocity increment AV r on the kth burn segment is given by 

K 

av l = AV i+l- AV i (3 - 1 ' 5) 


The burn time on the kth segment is then given by 


= T, - T. 
b i+i b i 


m c 
o 


AV. 

i 


AV,. 


i+1 


- e 


or. 



( 3 . 1 - 6 ) 
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(3.1-7) 


Defining the mass at time to be 


-AVj 


m. = m e 
i o 


the kth burn time can be written as 


m. 


-AV 


A 


T b t = TT c '^ 1 ' e C / 
k 


(3.1-8) 


It is now clear that, since can be expressed in terms of the 

equivalent velocity impulse on the kth segment AV , the optimal 

L k 

burn schedule, T, , T. , ... T, , can be obtained from the 
D 1 °2 °N 

velocity impulse schedule, AV T , AV, , ... AV , which 

1 i 2 TNT 

minimizes the characteristic velocity loss (3.1-2). 


The optimal N-hurn time-open burn schedule can be found by 
performing a parameter optimization on the impulse schedule 
and satisfying the appropriate constraints. For a specified V 
the sum of the N velocity impulses must equal the total required 
impulse AVj given by (2. 3-4). Furthermore, each burn must be 
"short", and the periapsis velocity after the next-to-last burn 
must be less than the escape velocity, as described in 2. 2. 
Therefore, the time-open multiorbit parameter optimization 
problem for specified N and V^ can be expressed as follows: 


Minimize: 0 


N 

C »I A \ f kV 2 

k=l * K 


(3.1-9) 
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Subject to: 


N 


1) Y AV t = AV t = \lv 2 ~ V 2 - V (3.1-10) 

i-i I, I 05 O 


k=l 


2 ) T < T b = V 

k max 


R l 




(3.1-31) 


3) V N-1 < V E 


(3.1-12) 


where 0 defines the cost function AV and V^__^ corresponds to 
the absolute periapsis velocity after the N-lst burn. 


The number of burns is not taken as a parameter in the above 
formulation, since 0 is not differentiable with respect to N. It 
can be seen that the cost function for the time-open maneuver is 
a minimum for a burn schedule consisting of N infinitely small 
impulses, which correspond to an infinitely large number of burns. 
However, the short burn and velocity constraints put an upper 
limit on the number of burns for a given V . Nevertheless, if 
V w is small — > 0 ), N will get very large, and the solution is 

clearly impractical. Therefore, the optimal time-open solution 
would have to be chosen on the basis of a transfer-time/perfor- 
mance-penalty tradeoff between solutions corresponding to various 
values for N. 


An approximate solution for the optimal time -open burn 
schedule can be obtained by assuming that the mass decrease 
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over each burn is negligible (e. g. , for high exhaust velocity 
powerplants). The burn time for each segment is then given by 

- e c J (3.1-13) . 

As a simple case, consider the two-burn transfer for which the 
cost function can be written as 



“2 ’ C . ( AV I, f l \ 2 + f 2 ) < 3 - M4 > 

Assuming that the inequality constraints (3.1-11, 12) are satisfied, 
the minimization of 0 2 i- s subject only to the velocity constraint 
(3. 1-10) given by 




AV, 


- (* v i. 


+ AVj ) - 

1 V 


(3.1-15) 


The solution to this simple two-parameter minimization problem 
can be obtained by defining the Lagrange multiplier v as a penalty 
on the constraint and by solving the Euler-Lagrange equations 
which describe the necessary condition for a stationary point, 
i. e. , 


• ^AV V ' * t>) = 0 (i = 1. 2 , ...N) 


(3.1-16) 


L i 


Assuming that f^ = 1 for all values of k (which is a very good 
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approximation), the required partials of (3.1-16) can be computed 



(3.1-17) 


Equating the expressions for WC given by (3.1-17), the obvious 
result is the relation 


AV t = AV 
L 1 X 2 


(3.1-18) 



(3.1-19) 


hi other words, the approximate two -burn, time -open optimal 
burn schedule is characterized by burns of equal magnitude and 
duration. 


Since the moment correction factor f. on the kth burn is. a 

k 

function of the kth velocity impulse AV T , removing the assumption 

k 

that f^ = 1 does not change the above result. Moreover, this 
result can be extended to the general N-burn time-open transfer. 
With the same initial assumptions used for the two-burn maneuver. 
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the N-burn optimization problem can be written as follows: 


min: 0 


N 

l A \ p k T b „ 2 

k-i k k 


Subject to: 

* = 



0 


(3.1-20) 


The necessary conditions for a stationary point are given by 


3 0 v 3 $ 

JSVT 3AV 
i . 


0 



y 3 tl) 

dSVf 


(3.1-21) 


3 0 
32\V t 


+ 


y 3 i/) 
3AV 

a N 


0 


where, without further manipulation, the solution is clearly given 
by equal velocity impulses and equal burn times: 



AV 
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N 


31 


k = 0, 1, 2, . . . 


(3.1-22) 



The general result can be stated as follows: For an N-burn, 
time-open multiorbit transfer for which the mass decrease over 
each thrust phase is negligible, the optimum burn schedule is 
prescribed by specifying each perigee burn to be of equal 
duration. This result is conditioned, of course, on the satis- 
faction of the inequality constraints. Hence, for a given V , 
a lower limit on the number of burns is established by the short- 
burn constraint (3.1-11), and an upper limit is placed on N by the 
escape velocity constraint (3.1-12). 

These results are summarized for a typical injection maneuver 

in Figures 3.1 and 3. 2. The approximate time-open simulation 

is described in Appendix C.l; the orbital parameters and power 

plant data used in all the simulations are listed in Appendix B. 

The vehicle selected for this study is a nuclear solid-core rocket 

with a capability consistent with the technology for the late 1970's 

or early 1980 's; it has been assigned a thrust to weight ratio of 

0.1 and a specific impulse of 800 seconds. Figure 3.1 shows 

the characteristic velocity loss AV ' as a function of the number 

of burns N for four different values of V . The endpoints for 

each family of solutions are specified by the short burn and 

escape velocity constraints as described previously. Figure 3.2 

illustrates the transfer time T as a function of N for the same 

s 

values of V b . It is interesting to note that, as the upper limit 
on N is approached for a given , the time in orbit asymptot- 
ically approaches infinity. The limiting case for each value of 
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V ro corresponds to an (N-l)-burn parabolic escape for which the 
time in orbit is given by the lower family of solutions corre- 
sponding to V ro = 0. The limits on N for each value of are 
listed in Table 3.1. 


Table 3.1 

Limits on N for the Approximate 
Time -Open Optimal Transfer 
(F/W = 0.1, I = 800 sec.) 


Lower Limit Upper Limit 


V » (fps) 

T > T 
b, b max 

k 

V N-1 =* V E 

0 

4 

CO 

5,000 

4 

31 

10, 000 

5 

8 

15, 000 

- 

- 


An upper limit on is specified by the short burn constraint. 
This fact accounts for the null solution corresponding to 

= 15,000 fps in Table 3.1. For the approximate time-open 
solution, the upper limit on the excess hyperbolic velocity is 
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given by 


V = VV„ - V‘ 

co max H E 

p max 


+ AV t - V‘ 
E I max/ E 


(3.1-23) 


where,, from (3.1-8), the maximum allowable velocity impulse 

AV t is given by 
I max s J 


AV. 


I max 


c In ( 1 - 


TV F 
b max 

m c 
o , 


= - c £n tl - 


F . B 


m c ' Li 
o ^ 


(3.1-24) 


For the maneuver simulated, these values are as follows: 


T, = 83 7 sec. 

bmax 


AVj = 2, 846 fps. 

xmax 


V = 14, 640 fps. 

<» max r 


These results, as well as those obtained for the limits on N, 
imply that the upper family of solutions shown in Figures 3.1 and 
3. 2 are invalid, and demonstrate that the short burn constraint 
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is somewhat restrictive in terms of final obtainable energy. 

Since the mass-loss approximation used in the above analysis 
rarely applies to typical propulsion systems, the approximate 
time-open solution is of little practical value. Nevertheless, the 
results are meaningful in the sense that they provide a simple 
heuristic framework from which to analyze the effects of and the 
restrictions imposed by velocity and burn time constraints on the 
optimal N-burn multiorbit maneuver. 

3. 2 Optimal Time-Open/ Time- Fixed Solution 

Consider the N-burn multiorbit injection maneuver for which 
the total transfer time is specified, that is, the time between 
the initiation of the first burn, which takes the vehicle out of 
parking orbit, and the escape burn, which occurs at the completion 
of the last intermediate coasting orbit, it assigned a particular 
value. Considering only the time constraint, the minimization 
problem can be stated as follows : 

N 

min 0 = C Y AV t f, T, 2 (3. 2-1) 

k=l K K 

subject to: 

*1 ” T s - T sd ■ 0 < s - 2 ' 2 > 

where T gd is the specified, or desired, transfer time and T s is 
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the actual transfer time, which is equal to the sum of the periods 
of the (N-T) intermediate orbits. The above problem can be 
restated as an unconstrained problem by defining the augmented 
cost function J such that 


J = 0 + (3. 2-3) 

where v, is a Lagrange multiplier which assigns a penalty to the 

unconstrained cost 0 proportional to the deviation of the actual 

time T from the desired time T Hence, if v, = 0, the 
s sd 1 

augmented cost J equals the unconstrained cost 0, and the 
minimization corresponds to the time-open maneuver. 

The optimal burn schedule for the time-fixed maneuver could 
be obtained by means of a conventional constrained parameter 
optimization algorithm, such as a gradient projection scheme. 

An alternate, simpler approach will be considered here. The 
major advantage of the method to be derived is a reduction in 
the dimension of the problem to a form which lends itself to 
simple iterative methods. For purposes of simplification, the 
inequality constraints described in the previous section will be 
temporarily ignored. 

First, let the variable AV . be the total velocity impulse 
effected up to the total burn time T^ , which for constant 
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thrusting is given by 



(3.2-4) 


Recall the segment variables for the kth thrust phase previously 
defined as 



AV i+1 - AVj 



(3.2-5) 


where k varies from 1 to N. The problem can now be restated in 
terms of the variables AVt and , which behave as continuous 
variables over the total burn interval of a multiburn-multiorbit 
trajectory. The characteristic velocity loss can be rewritten as 


N-l 


AV * = c . I ( AV i+ i - AV i) ( T b. , - T b. ) (S - 2 - 6 > 

l+l l 


1=0 


where f ^ is a function of the quantity (^AV^ - AV^. The total 
time in orbit is related to the variable AV^ by 

N-l N-l 


T s = l T i ( AV i) = l [ 


R 


n 3/2 


i=l 


i=l 


2(1 - H (v o + AV t ) ; 


(3.2-7) 


where T^ is the period of the orbit following the ith burn. 
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The necessary condition for a stationary point of the augmented 


cost function J is given by 


a j 

asv: 


a (0 + v 1 0^ 


(3.2-8) 


Assuming that f ^ = 1 and taking the required partials of AV 
given by (3. 2-6), the necessary condition becomes 


-T*. + T 


V +2 AV t T, ^ + 2 AV t T, v - r . w •■■ ■ 
b k-l \ b k d ^ V i l k-l b k-l 5 ^ v i J 


( AV i) . 


l a 


(3.2-9) 


where the k subscripts denote the segment quantities defined by 
(3.2-5). Utilizing (3.1-6), the required partials of the segment 
burn times are found to be 


aAv. 


-AV. 

c 


-AVi 

m 

o Q c 

~F~ e 


If the quantity is defined, where 

( 

Q, = \l-e c / 



(3.2-9) can be written as' 



+ 



/m .\ 2 

2 \ t ) «k A \ 



«k-l AV I 


k-1 



aT s ( AV i) 

bAvr- 


= o 


(3.2-12) . 


Rearranging, (3.2-12) can be reduced fco 




+ 2AV, 


+ 





k-1 


+ 





+ 


C c 

CO 


a T s ( AV i) 

§"AV i 


0 

(3.2-13) 


The variables AV.. corresponding to the optimal burn schedule 
must satisfy (3. 2-13). The optimal burn schedule could be 
obtained by simultaneously satisfying (N-l) equations of the form 
(3. 2-13) and the velocity constraint on the last burn given by 


or. 


V + AV = V - VV 2 + V 2 
o A N H p v E 


AV N = VV E 2 + V » 2 - V o = AV I (3 - 2 ' 14) 


A simpler method for obtaining the optimal N-burn thrust 

program can be arrived at as follows: First, expand the kth 

41 



burn time given by (3.1-6) in a Taylor series about A-V^. To 
first order, this expansion is given by 


vr v * +A \ r 1 n +A v r> 


or. 


T. 


b i+i 


T b t + (^1+1 - Av y 


\ m 


-AV. 


(3.2-15) 


Hence, the duration of the kth burn is approximated to first order 

by 


m. / 

T, - t, « AV.,, - AV. 

b i+l b i F \ 1 1 *7 


which, in terms of segment quantities, is given by 



(3. 2-16) 


Substituting this approximation for the burn time into (3. 2-9), 
the first order necessary condition reduces to 



(3. 2-17) 
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Recalling the definition of AV T given by (3. 2-5), (3. 2-18) can be 

A k 

utilized to give an expression for the total velocity impulse A\h + , 
at the end of the kth burn as follows: 



(3.2-19) 


The partial derivative of the transfer time, given by (3. 2-7), can 
be expanded as follows: 



Equation (3. 2-19) gives the useful result that the total velocity 
impulse AV. + j at the end of the kth burn can be calculated from 
a knowledge of the velocity impulse AV T on the previous 

Vl 

segment, the total velocity impulse AV. at the beginning of the 
burn, and the value of the Lagrange multiplier Vy In other words, 
given the initial velocity impulse and Vy an optimal trajectory is 
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completely specified by (3. 2-19). Furthermore, the validity of 
the use of a first-order approximation for the burn time is 
ensured by the fact that it is contained in a first order necessary 
condition. 

The N- dimensional parameter problem has thus been reduced 
to one of two independent variables, AV^ and u. . By means of an 
iterative procedure, AV^ and v-> can be chosen such that a unique 
transfer time can be obtained in exactly N-burns. An optimal 
•N-burn schedule results if the burn time and escape velocity 
constraints are also satisfied. A desired final energy, i. e., 
a specified V ro , can be achieved for a fixed-time transfer by 
iterating on AV^ until the total impulse on the Nth burn is such 
that (3, 2-14) is satisfied. 


The recursion formula (3. 2-19) can likewise be used to 
generate optimal time- open burn schedules by setting the 
Lagrange multiplier to zero in (3. 2-19), i. e. , 


Av i+1 


AV. + AVj 


k-1 




/ra i~l 
\ m i / 


^ 2 


+ 2 


v m iy 


(3.2-21) 


This time-open recursion can also be expressed as follows: 





(3. 2-22) 
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in other words, for the optimal time-open transfer, the velocity- 
impulses on succeeding burn segments of the equivalent impulsive 
trajectory are related by a function of the mass ratio between 
the two segments given by 


m. *i 
i-I 

m- 


A y i - A y i _ 1 


-av t 


k-1 


= e 


(3. 2-23) 


Therefore, for constant thrusting over the burn segment, the 
equivalent velocity impulses on neighboring segments of a 
time-open maneuver are' : r elated by a function of the velocity 
impulse on the earlier burn. 


Equation (3. 2-22) gives a measure of the error in the 
approximate time-open solution discussed in 3.1. It verifies that, 
for a negligible mass decrease over each burn, the ratio given 
by (3. 2-23) approaches unity, and neighboring burn segments have 
equal burn times. For a specific impulse of 800 seconds 
(c « 2.6' 10^ sec. ), the velocity impulses would typically have 
to be less than 300 fps (which gives a value of 0. 99 for (3. 2-23)) 
to yield a valid approximation. Clearly, such small velocity 
increments would require large values for N and total transfer 
times out of the range of practical concern. 

The following chapter will describe a technique that can be 
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readily employed to predict optimal time-open and time-fixed 
multiorbit transfers using the recursive formulation derived 
above. 
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CHAPTER IV 


ANALYSIS OF A TYPICAL INJECTION MANEUVER 


4.1 Time-Open Solution and Determination of Initial Values 

The optimal N-burn thrust program for specified values of V ro 

and T can be found by iterating on the parameters AV. and v , , 
s 11 

where each intermediate impulse schedule AV T , AV T , ... j\V , 
is generated by the recursion formula (3. 2-19). For the time-open 
case (v^ = 0), the optimal N-burn schedule (for a a specified VJ 
can be computed by iterating on AV^ until the sequence generated 
by (3. 2-21) yields a value for AV^ that satisfies the final 
velocity constraint given by (3.2-14). The impulse schedule 
obtained corresponds to the optimal time-open N-burn thrust 
program if the short-burn constraint (3.1-11) is satisfied and if 
the solution is constrained to exactly N burns by the escape 
velocity constraint (3.1-12). 

' If the appropriate value for the initial velocity impulse AV-^ is 
selected, a solution corresponding to the desired value for N can 
be obtained. However, if a "good” initial value is not selected, 
a burn schedule corresponding to a different N will be generated 
and Newton steps will not cause the sequence to converge to the 
desired N-burn solution. In other words, solutions corresponding 
to different values of N are independent and self-contained. This 


47 



is evident from the nature of file velocity constraint (3.1-12); each 
N-burn solution is separated from its adjacent (N-l)-burn 
solution by an (N-l)-burn parabolic escape solution. The 
independence of N-burn solutions is illustrated in Figures 4. 1 
and 4. 2 , which are the results of the simulation described in 
Appendix C.2. These figures plot, respectively, V and T 
against the parameter AV^ for the time -open, velocity unconstrained 
injection maneuver specified by the parameters listed in Appendix 
B. Not only do these plots demonstrate the independence of N- 
burn families, but they can be utilized to give excellent starting 
values for AV^ for a particular N-burn velocity- constrained 
solution. In other words, if the initial impulse A.V^ corresponding 
to a specified V and a particular value for N is selected from 
Figure 4.1, it is guaranteed that (3,2-19) will generate an impulse 
sequence very close to the desired one, and the Newton iteration 
will converge in a few steps to the optimal N-burn schedule. 

Before discussing the effect of the time constraint (nonzero v^), 
consider the difference in the optimal time-open solution generated 
by the above method as compared to the approximate time-open 
solution derived in 3. 1. Figures 4. 3 and 4. 4 present plots of 

•JL* 

AV‘ vs. N and T vs. N, respectively, for the time- open 

D 

solution. The locus of solutions corresponding to the approximate 
time-open case is shown on each figure to provide a simple means 
of comparing the two methods. The difference in the magnitudes 
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Figure 4.1 V m vs AV^ For Time-Open, Velocity-Unconstrained Maneuver 







Figure 4. 2 T vs AV. For Time-Open, Velocity-Unconstrained Maneuver 

S X 




T 

s 

(hours) 


i 

I 


I 

i 

i 

i 



Figure 4.4 T vs N For Time-Open Maneuver 

S 
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between the values given by each method is not surprising, as 

the exhaust velocity for the powerplant used in the simulation 

( 

is not particularly large. Nevertheless, these results do show 
that the recursive solution gives significantly better results for 
an injection maneuver like the one simulated. Furthermore, 
contrary to the approximate results given in 3.1, there exists an 
optimal time-open solution for a final velocity of 15, 000 fps 
corresponding to a 5-burn multiorbit trajectory. Therefore, 
limits on N for specified values of and on for specified 
N's can be obtained from the results of the time-open, velocity- 
unconstrained simulation. Although solutions only for values of 
N less than 10 are illustrated in Figures 4. 1 - 4.4, the same 
arguments apply for large values of N. This study will ignore 
large JSF solutions, however, as the transfer times and trajectories 
associated with them are too large for practical purposes. 

4. 2 Effect of the Time Constraint on the Optimal Burn Schedule 

Recalling the necessary condition for a stationary point given 
by (3. 2-8), the Lagrange multiplier v, which can be written as 

- B0 /SAV. 

V = wjr (4.2-1) 

represents the sensitivity of changes in cost with respect to changes 

in the value of the constraint; or, in a plot of 0 vs. ip, v represents 
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the negative slope of the curve. If the stationary point is a local 
minimum, values of >p different from zero must be assigned a 
positive cost. Hence, v must be positive for positive values of 
ip and negative for negative ip. This relation is sketched in 
Figure 4, 5 



Figure 4.5 0 vs Ip 

Applying this argument to the multiorbit- injection problem, it is 
clear from (3. 2-2) that positive corresponds to T g > T g ^ and 
negative v ^ to T g < T Thus, a sketch of the characteristic 
velocity vs. the time in orbit for a specified and a particular 
value for N should resemble the curve shown in Figure 4.6, where 
the minimum value of AV' corresponds to the optimal time-open 
performance penalty and corresponding transfer time. 
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A physical interpretation of the effect of the time constraint 

could be given as follows: For the time-open case = 0), 

each succeeding burn segment will increase in length by an 

amount proportional to the velocity- increment ratio given by 

(3.2-22). Since nonzero values of v affect the lengths of the AV T 

increments on succeeding segments through the second term under 

the radical in (3. 2-1 9), a positive value for v ^ tends to make the 

AV t increase or successive burn-segments greater; whereas, a 

negative reduces the rate of impulse expansion on successive 

burn segments. In other words, for v, > 0, (N-l) AV T segments 

k 

are "stretched-out" to begin and end at higher absolute periapsis 
velocities, which correspond to intermediate orbits of greater 
periods and to a longer overall maneuver time T . Similarly, 
for pj < 0, burn begin times are scheduled at periapse velocities 
less than those of the corresponding time-open maneuver, and the 
total injection time T g is less than that for the time- open case. 
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The effects of a time constraint on a typical maneuver 

(Appendix B) are illustrated by Figures 4. 7 and 4. 8, which 

summarize the results of the simulation described in Appendix 

C, 3. Figure 4, 7 is a plot of the specified transfer time T vs. 

s 

the characteristic velocity loss AV for a 5 -burn injection to a 
specified final velocity of 5000 fps. Note that the Lagrange 
multiplier used in the simulation, X, is weighted by the constant 
C m and takes on a sign opposite to that shown in Figure 4. 6 (since 
the constraint is defined with the opposite sign in the simulation). 
In other words, the Lagrange multiplier' shown on Figures 4. 6 
and 4. 7 is given, for convenience, by 


X 


(4.2-2) 


Each point on the curve in Figure 4. 7 corresponds to an 
optimal time-fixed, velocity-constrained burn schedule generated 
by the recursion formula (3. 2-19). Each solution was generated 
by selecting a value for X and performing a Newton iteration on 
the variable AV^ until the final velocity converged to within 5 fps 
of the desired value. The multiplier X was varied in the positive 
and negative direction from X = 0. AV^ was initially selected 
for the time-open case with the aid of Figure 4.1, and each 
consecutive point on the curve was determined using a greater 
value for X and the value of AV^ for which the solution converged 
on the previous point, TheX recursion was terminated for positive 


57 



values when the lower limit on T g was reached, i. e. when the 
short burn constraint was exceeded. The recursion on negative 
X was terminated when AV^ increased to a value such that the 
optimal recursion (3.2-19) could not be constrained to 5 burns. 

The endpoints were not determined precisely, as they are not 
required in terms of desired results. 

The predicted effect of the time constraint on the schedule of 

velocity impulses is also confirmed by the simulation. The 

stretching and compressing effect of X on the impulse sequence is 

illustrated in Figure 4. 8. The Figure compares the thrust 

program for the optimal time-open 5-burn injection with two 

corresponding optimal time-fixed transfers, one with a transfer 

time greater than the time-open maneuver and one with a shorter 

transfer time. For each -of the three cases given, the total . 

required impulse AV^ (for an injection to V OT = 5000 fps) is 

divided graphically into 5 burn segments. In each segment the 

equivalent velocity impulse AV, and the burn time T, are 

k D k 

given in units of feet-per- second (fps) and seconds (sec), 
respectively. The velocity impulse variable AV^ can be read 
from the scale at the bottom of the figure. The semimajor axis 
(in miles) and the period (in hours) of each intermediate coasting 
orbit are listed at the end of the appropriate burn segment for each 
case. The flight path profile for the optimal time-open 5-burn 
injection corresponds to that shown in Figure 1. 
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Figure 4. 8 Effect of Time Constraint on Burn Schedule 

















4.3 Determination of the Optimal Time-Fixed Burn Schedule 

For a multiorbit injection maneuver of specified final energy, 
a tradeoff must be made between the transfer time and the 
performance penalty. Moreover, in order to obtain the optimal 
burn schedule, maneuvers corresponding to different values of 
N must be compared. This tradeoff can be readily accomplished 

■T* 

with the aid of a plot of T vs &V for various N-burn solutions 

s 

satisfying a common final energy requirement. For the injection 
maneuver considered previously, this type of analysis can be 
carried out with the aid of Figures 4. 9 and 4. 10, which plot 
T g vs AV ' for N-burn families corresponding to V a = 5000 fps 
and V m = 10, 000 fps, respectively. Each N-burn curve is gen- 
erated as described for the 5-burn solution given in Figure 4. 7; 
however, for large values of N, the lower limit on each curve 
is specified by their intersection with the (N-l)-burn curve rather 
than by the short-burn constraint. The lower limit on N for a 
given family of curves is established by the short-burn constraint. 
N is not taken to its upper limit as these solutions are not of 
interest. 

It is interesting to note the double solutions given by the inter- 
section of curves for larger values of N, The values of T g and 
AV corresponding to one of these points represents two different 
N-burn optimal transfers; hence, some other criterion must be 
utilized to select the better solution. In general, the optimal 
N-burn solution for an injection maneuver of specified should 
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lie on the lower edge of the profile described by the family of 
curves. 

If a particular transfer time is specified, the simulation used 

to generate the T vs. AV" curves will not yield the optimal 
s 

burn schedule in one convergence cycle, since the value of 

that corresponds to a specified T gc j is not known. If a particular 

value of T is a critical mission requirement, the optimal N-burn 
s 

schedule can be obtained by a two-dimensional iteration on the 

variables AV^ and \. A two-dimensional Newton iteration which 

yields a particular time-fixed, velocity-constrained, optimal 

N-burn schedule is described in Appendix C.4. This technique 

has several disadvantages, however. First, it requires initial 

values of AV^ and X which guarantee that the solution will locally 

converge to the desired N-burn maneuver. Second, the solution 

obtained is not a "global" one; hence, it does not guarantee the 

best thrust program. Therefore, the two-dimensional search 

method should be used as a secondary routine to obtain a particular 

fixed-time burn schedule in a local N-burn region. In other words, 

the time/cost tradeoff should be made from the T vs. AV' curves, 

s 

then the optimal time-fixed solution can be generated by a two- 
dimensional search routine which utilizes initial values obtained 
from the one dimensional simulation. 


The one- dimensional simulation is not without its difficulties, 

however. For example, the optimal recursion will not converge 
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for a parabolic escape maneuver. This tendency can be explained 
by referring to Figure 4. 1, where it can be seen that the slope 
of each N-burn curve is infinite at = 0. This implies that 
the initial value for AV^ must be exact for the solution to 
converge. Also, for large values of and N, the Newton step 
size must be heavily controlled to constrain the solution to N 
burns. Fortunately, however, the latter concern is outside the 
range of interest of this analysis. 

4. 4 Accuracy of the Solution 

In order for the characteristic velocity loss equation to be 
valid, the burn time must be sufficiently short and the variables 
r and 0 must vary very little during each thrusting phase. A 
measure of the accuracy of the solution obtained in this analysis 
is provided by the values taken on by some of the thrust parameters 
which behave as error quantities. The most obvious measure of 
the validity of the solution is the burn time. Longer burn times 
give a less accurate value for the performance penalty, and burn 
times greater than the upper limit given by {3. 1-11) deem the 
solution invalid as far as this analysis is concerned. Another 
measure of accuracy is the magnitude of the radial displacement 
effected during each burn {equation 2.1-4). This displacement, 
which must be negligible with respect to the radius R of the 
initial orbit, provides an indirect measure of the span of the 
thrust arc. These two quantities are tabulated below for various 
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values of N corresponding to the optimal time-open burn schedules. 
Also tabulated i's the moment correction factor, which was 
assumed to be unity in the derivation of the recursion formula 
(3.2-19). For each case, the value of the error quantity repre- 
sents the "worst case" value for that particular N-burn solution. 
For the time-open maneuver, the worst case, i. e. , the largest, 
values of TV and d occur on the initial burn segment, whereas, ■ 
the worst case value for f^_ i. e. , the one that deviates the most 
from unity, occurs on the escape segment. 



V = 10, 000 4 ■ 788 22.86 i . 99972 

5 ! 633 | 11. 52 j .99981 

6 I 528 I 6.60 | . 99987 

7 453 j 4.12 ! .99990 

8 397 ' 2.75 j .99992 
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For the N-burn time-open maneuver considered in Table 4.1, 
the burn times stay well within the short-burn- constraint but 
diverge toward T. (837 sec) for small values of N, part- 

icularly for larger V . The radial displacement at perigee 
likewise gets worse for low values of N and large V^; however, 
even for the worst N-burn case shown, the magnitude of the 
displacement is only 0.56% of the radius of the initial orbit. 

The values of f^, listed are all within at least 0. 01% of unity; 
thus, the approximation used in the derivation is confirmed by 
the results of the simulation. 

It is also useful to note the variation of these error quantities 
along the segment of an optimal N-burn multiorbit transfer. For 
the three cases described in Figure 4. 8, the error quantities 
propagate as shown in Table 4. 2. It is interesting to note that 
for the time-open case (X = 0), the error quantities vary little 
over the trajectory. In fact, the displacement is approximately 
the same at each perigee. In a sense, one could conclude that 
these quantities provide a relative measure of the characteristic 
velocity loss along the trajectory. The only trend indicated by 
Table 4. 2 that is of any concern is the rapid divergence of the 
error quantities on the last segment of time-fixed maneuver 
corresponding to positive X. The values listed on Table 4. 2 for 
this case are within the limits for a good approximation; however, 
it should be realized that, for some multiorbit injection problems. 
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Table 4. 2 

Variation of Error Quantities along 
as Optimal 5-Burn Multiorbit Trajectory 


(F/W = 0.1, 

Segment 

Number 

k 

I = 800 
s 

P 

Burn 

Time 

T. 

b k 

(sec) 

sec, V ra = 5000 fps) 

Radial Correction 

Displacement Factor 

-k f k 

(miles) 

-5000 1 

610 

10. 55 

. 99990 

2 

592 

10. 45 

. 99989 

3 

572 

10.18 

. 99987 

4 

546 

8.42 

. 99986 

5 

454 

5.12 

.99987 

0 1 

587 

9. 38 

. 99990 

2 

5 71 

9. 37 

. 99989 

3 

555 

9. 37 

. 99988 

4 

538 

9. 37 

. 99987 

5 

521 

9. 36 

. 99985 

20, 000 1 

542 

7. 43 

.99992 

2 

534 

7. 74 

. 99991 

3 

531 

8.44 

. 99989 

4 

.540 

11.19 

. 99987 

5 

626 

18.12 

. 99979 


the approximation may degenerate for large values of and 
low values of N. 


Other errors associated with the analysis described in this 
study are the numerical errors inherent in the Newton iteration. 
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However, for the one- dimensional search utilized, these errors 
can be made negligible by specifying a more restrictive 
convergence criterion on the error associated with the final 
velocity constraint. The analysis has demonstrated that the 
Newton iteration will converge in a few steps for specified 
tolerances of less than 0. 5 fps on V ro , i. e. , errors less than 
0. 01% of the magnitude of V . Hence, numerical errors can be 
made small enough to neglect. 

Finally, it should be recalled that there may be a small error 
in the solution, due to the use of a first-order approximation for 
the burn time. This type of error can be avoided by finding the 
impulse schedule that exactly satisfies the first order necessary 
condition given by (3.2-13). Appendix C. 5 outlines an N-dimensional 
search routine that is designed to do this. The results of this 
simulation are not presented, as the solution could not be driven 
to within tolerances that justify its use as a measure of validity 
of the first-order burn time approximation. Furthermore, an 
exact solution of the necessary condition may not be justified in 
light of the nature of the approximation for the characteristic 
velocity loss. 
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CHAPTER 5 


CONCLUSIONS 


This study has presented the followsing results: 1) A simple, 
closed-form method for calculating gravity losses on multiburn- 
multiorbit finite thrust trajectories can be developed utilizing 
Robbins’ "impulsive approximation". 2) This approximation can ‘ 
be applied to the multiorbit injection problem for which a simple 
algorithm can be derived to recursively generate optimal N-burn 
thrust programs. 3) The recursion algorithm can be used to plot 
families of time/cost curves for a specified value of the excess 
hyperbolic velocity, and these curves can be utilized to select 
initial values from which the optimal N-burn time-fixed injection 
maneuver' can be generated. 

The numerical values presented in this study correspond to a 
particular injection maneuver. It is important to note that the 
results of the simulations are very sensitive to the selected 
values of the specific impulse and thrust-to-weight ratio. Clearly, 
higher values for P/W could be utilized, and the equivalent 
injection maneuver could be accomplished in fewer burns or with 
a smaller penalty, ha any case, however, due to the nature of the 
short-burn criterion, the thrust parameters selected will impose 
restrictions on the final obtainable velocity and the minimum 
allowable number of burns, ha spite of these restrictions, the 
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method presented in this paper provides a simple and efficient 
means for predicting optimal thrust programs for N-burn 
mulfciorbit injection maneuvers and provides a framework from 
which optimal burn schedules can be derived for other multiorbit 
missions. 
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APPENDIX A 


DERIVATION ■OF CHARACTERISTIC VELOCITY LOSS 
OVER A SEGMENTED FINITE- THRUST TRAJECTORY* 

The differential equations of motion for a vehicle in a vacuum 
are given by 

r = v 

v = g (r, t) + a(t) (A-l) 

where r and v are, respectively, the position and velocity vectors 

describing the vehicles motion, g. is the gravity vector, and a 

is the thrust acceleration vector. If r and v are perturbed about 

2 

a nominal trajectory, such that terms of order 6r and higher 
can be neglected, A-l can be expressed by 

6r = 5v 

6v = G6r + 6a {A- 2) 

where G is the gravity gradient matrix. If there is no thrusting, 
i. e. , if the vehicle is coasting under the influence of gravity, the 

* This derivation closely follows that of Robbins [12]. 
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perturbed equations of motion can be written as 


6r1 


O 

I 


"6r ” 

6v 


G 

o_ 


8v_ 


or, in second order form, 

6r = G 6r 


(A- 3) 


(A-4) 


Consider the adjoint system to (A-3) given by 


'if 


0 

-g" 


V 

_X 


-I 

o_ 


-A- 


(A- 5) 


where £ and X are (3 x 1) vectors adjoint to £ and v, respectively. 
In second order form, (A-5) becomes 


X = GX (A- 6) 


Thus, X, the adjoint to the velocity vector or "primer vector" of 
Lawden, satisfies the same differential equation as the perturbed 
position vector. 


It can be shown that the primer vector satisfies the identity 

* 

X * 6v - X * 6r = constant (A- 7) 
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For non- zero thrust over the time interval (t , t f ), it can 
similarly be shown that 

[x * 6v - X • 6r]* f = X • 5a dt (A- 8) 

o J o 

recognizing that (A- 7) takes on a different constant value outisde 
the burn interval. The variation in the characteristic velocity 
from its nominal value is given by 

* *fc "t 

6(AV) = fx ' 6v - X ' 6rl, f + { f (6a" - X * 5a) dt 

L “ ~ - ~o JT o ” " (A- 9) 

where 6a' is the variational ’thrust acceleration over the nominal 
finite-thrust interval. Assuming that both trajectories start 
from the same initial state, i. e. , 6v (t ) = 6r (t Q ) = 0, and that 
they both satisfy the same final conditions (such that 6r (t^) and 
6v (tj.) are small enough to justify linearity), X (t^) and X (tp 
can be chosen such that (A- 7) is identical to zero over the 
interval. (A- 9) then reduces to 

6(AV) = f. f (6a' - X • 6a) dt (A-10) 

Jt o 

If the nominal trajectory is optimal, the right hand side of 
(A-10) must be non negative for all perturbations of the acceler- 
ation. If optimal steering if utilized, that is, if the acceleration 
is applied in the direction of the primer vector, X » _6a = j X | 5a', 
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and (A-10) becomes 


6(AV) 



(1 - jxl ) 6a" dt 


(A-ll) 


which implies that a necessary condition for an optimal impulsive 
trajectory is that, at the time the impulse is applied, the 
magnitude of the primer vector must be unity, and for all other 
t in the interval, j xi <1. 


Let AVj be the optimal impulsive velocity on a time uncon- 
strained trajectory, that is, a trajectory on which the time of 
occurence of the impulse is determined by local optimality with 
respect to fuel consumption {with the time constraints playing no 
role except possibly to exclude other local optima). Also, let 
AVp be the characteristic velocity over a finite thrust trajectory 
which is not necessarily optimal, but which is close to the 
optimal impulsive trajectory with respect to position deviations 
and which satisfies the same end conditions. Then, if is 
the thrust acceleration vector on the finite thrust trajectory and 
a is the angle between and X, (A-10) can be written as 

c tf 

AV f - AVj = ^ (a F ' - X ' a p ) dt 
t o 


(1 - Ixl cos a) a„ dt (A-12) 

t ~ 1 

o 
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If are the times of occurence of impulses over a multiburn, 
segmented trajectory, i. e. , one that can be broken up into a 
series of short burns about local optima and separated by coasting 
arcs over which a^ = 0, the integral (A-12) can be expressed as 
a sum of integrals, each taken -over a short burn period near its 
respective Furthermore, cos a. can be expanded as 

cos a. « 1 - vjr o: 2 (A-13) 


and the primer vector magnitude can be expanded in a second 
order Taylor series expansion about t^, i. e. , 


\(t)| = |x(t k )| + |i(t k )| <t-t k ) +\ |x(t k )| (t-t k ) 


k 2 


= 1 + \ “-‘k ’ 2 


(A-14) 


since t k is at a local maximum. Furthermore, since |X(t k )| < 0 
at the local maximum, (A-12) can be expressed as 

n ^ek 

AV F - AV l- l i [I “ 2 - I < k -V 2 lik >F dt 

k=l t bk (A-15) 


which is valid for a segmented trajectory of n short burns, each 
of which extend between their respective burn begin times 
and end times t gk . 
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If t , is the centroid time of the thrust acceleration profile 
on the kth burn, the first moment of the thrust acceleration 
about its centroid is defined to be zero, i. e. , 

^ek 

Mjk = $ (t-t c k> a F tt) dt = 0 (A-16) 

*bk 

and the second moment is given by 

— e k 

M 2k C \ a-t ck )2 a F (t) dt (A ' 17) 

" 'bk 

Assuming that a varies approximately linearly during a burn, i. e. , 

a(t) ~ a ck + (t-tfc) a C k (A-18) 

and if the burn time is referenced to the centroid time, i. e. , 


t-t k (t ""W + ( 


(A-19) 


the characteristic velocity loss (A-15) can be written as 


n 


AV P-* V I “ y 1 U“ok 2 -('ck - *k T\] AV t 


k=l 


+ 


L“ck 


lilkKJ* 


(A-20) 
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where the velocity increments AV^, > 

^ek 

AV k = V a F '<« dt 

T bk 


(A- 21) 


are given to sufficient accuracy by the impulses of the optimal 
impulsive trajectory. It can be seen from (A-20) that the optimum 
choice for the centroid time t ^ is the time of occurence of the 
impulse t. ; and, it can be- shown that this choice for t^ is that 
which also results by requiring the component of position 
deviation parallel to X (t^.) to be approximately zero both before 
and after the burn. Finally, with optimal .timing ,(t ^ = t^) and 

9 

optimal steering (a^ is parallel to X(t) such that a ^ = 0), 

(A-20) reduces to 


n 


AV 


F 


AV t 


k=l 


Ak 


M 


2k 


(A-22) 


where the second moment can be accurately approximated as a 

function of the velocity impulses AV T on the equivalent impulsive 

L k 

trajectory. By differentiating the identity 


X = 


(A-23) 


(A-6) can be utilized to obtain the relation 


A ' i + A ‘ G A = |X| j A i + |il 


(A-24) 
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Then, since |x| = 1 and |X| = 0 at t k , |X k J can be solved for as 

|Xj k - X • GX + jxj 2 (A-25) 

Substituting into (A-22), the characteristic velocity loss is given 

by 

n 

£V f - AVj « \ (-X * GX - |X| 2 ) M 2k (A-26) 

k=l 

Therefore, the characteristic velocity loss over a segmented 
trajectory consisting of n optimally steered and timed finite-burn 
segments of short duration can be computed in terms of the primer 
vector solution at the centroid of each burn and the velocity 
impulses required to effect the equivalent transfers. 
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APPENDIX B 


ORBITAL PARAMETERS AND 
POWER PLANT DATA 

The following data correspond to an injection maneuver from 
a circular orbit 100 miles above the surface of the earth using a 
nuclear solid-core rocket consistent with the technology of the 
late 1970's - early 1980's. These values are incorporated into 
all the simulations used in- this study. 

Orbital Parameters: 


M 

= 1.4076 

, n 16 p.3 , 2 

' 10 ft /sec 

R 

= 4.0632 

, n 3 

* 10 mi. 

V 

0 

= 2.561 

■ 10 4 fps 

T 

o 

= 5.262 

3 

• 10 sec. = 88 min 

V E 

= 3.622 

' 10 4 fps 


Powerplant Data: 

I g p = 800 sec. 

F/W =0.1 

c = 800 g = 2.57 * 10 4 fps 

Short-burn Constraint: 

T b = 837. 6 sec. 

D max 
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APPENDIX C 


COMPUTER SIMULATIONS ‘ 

To avoid notational difficulties and unnecessary repitition, a 
brief summary, rather than a complete program listing, of each 
simulation is provided in the following sections. The orbital and 
powerplant data used in each program correspond to the maneuver 
defined in Appendix B. 

C.l Approximate Time-Open Solution 

For a specified value of and N, the approximate time-open 
solution is a completely closed-formed calculation. The simulation 
repeats itself only for different values of N and V . The total 
required velocity impulse AV^ is calculated from (2. 3-4) and the 
equivalent impulse and the burn time are calculated from (3.1-22). 
For a given , the range of N is determined by the upper and 
lower limits prescribed by the escape velocity and burn time 
constraints given by (3.1-11) and (3.1-12), respectively. In 
addition to the segment quantities, this simulation computes the 
characteristic velocity loss using (3.1-2). The results of the 
simulation are summarized in Section 3.1. 

C. 2 Time-Open, Velocity- Unconstrained Simulation 

If is unspecified and is set equal to zero, the curves 
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given in Figures 4. 1 and 4. 2 can be generated by selecting 

various values of AV-^ between zero and the approximate time-open 

AV and by performing the recursion given by (3. 2-21) until 
unax 

the escape velocity is exceeded. ,The value for the index is set 
equal to N on this step, and V ro , is computed from {3. 2-14). Hence, 
the array AV^, AVg, . . . AV^ is defined, and all orbital quantities 
(m., a., T.) and burn segment quantities {AV , T, , d. , f, ) can 

111 * JLj^ K K 

be calculated using the equations given in Chapters 2 and 3. 

C. 3 Time-Fixed Recursive Solution for Specified V w 

If V ra and X are specified and if a good initial value for AV^ is 
selected (from Figure 4. 1), such that the value of N given by the 
recursion (3. 2-19) and the escape velocity constraint (3.1-12) 
corresponds to the desired value, a one- dimensional Newton 
iteration on AV^ can be effected such that the final velocity given 
by (3.1-2) will converge to the desired value in a few steps. The 
classical one -dimensional Newton step is defined as 

Ax = (C. 3-1) 

where f(x) is the value of the function that is to be driven to zero 
by proper choice of the variables x and f'(x) is the derivative of 
that function with respect’ to x. 
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In this simulation f(x) corresponds to the final velocity- 
constraint, given by 

f(x> = *(AVj> = V wd - V o 

and f (x) is given by 

a^AV^ -9v ra 
f (x) “ T5Vy = 

The derivative is computed numerically by a balanced difference, 
i. e. , 


(C. 3-2) 


{C. 3-3) 


dV V (AV, + 6) - V (AV, - 6) 

co co £ CO X 

dSv[ ~ 2~5 


(C. 3-4) 


where 6 is a small increment, specified in this simulation as 


6 = .01 AV 1 


Subsequent values of AV^ in the iteration are computed as 
follows ; 


6 < Av i>k« 


< AV i»k 


a * (AV l , k 


(C. 3-5) 
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The iteration is continued until the velocity constraint (C. 3-2) 
has converged to within some specified value of zero. The 
convergence criterion used in the simulation is given by 

$(AV X ) < 5.0 fps (C. 3-6) 

This criterion gives good results as well as convergence in two 
or three steps. To assure convergence the step size has to be 
controlled in the region of the solution. The control used in this 
simulation is simply a halving of the kth step size 6 (AV'-^)^ if 
the error, given by (C. 3-2), is not decreased on the kth step. 

When the iteration converges to the optimal impulse sequence, 

o> 

the orbital and segment quantities as well as T g and AV' are 

computed, as in the velocity unconstrained simulation. The iteration 

cycle corresponds to particular values of V , N, and A. and 

represents one cycle of the iteration on X described in Section 4. 2. 

When a sufficient number of cycles are run to plot a complete 

T vs AV’ curve, new initial values for AV.. are read in and another 
s 1 

N-burn family is' simulated. The results of this simulation are 
summarized in Figures 4. 9 and 4. 10 for two values of . All of 
the orbital and segment quantities given in Chapter 4 for the 5 -burn 
maneuver are computed by this simulation. 
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C.4 Two-Dimensional Search for Specified V and T 

If the transfer time T g is specified along with V m , the optimal 
N-burn schedule can be obtained by selecting the proper starting 
values and iterating consecutively on the variables AV^ and X 
until the velocity and time constants are satisfied. The iteration 
follows that described in C. 3; the two dimensional Newton step 
is defined below. 


Let the generalized 2-dimensional constraint $ have components 
^ and ^ where is the velocity constraint given by C, 3-2 and 
i [>2 is the time constraint given by 


= T sd - T s (C - 4 - 1) 

The step 6(AV- L ) J SX that drives (C.3-2) and C. 4-1) simultaneously 
to zero is given by the solution to the following equations: 

5 (AV i )+ -al 6(X) = -*1 

(C. 4-2) 

a ip 9 

d *2 = ^ 6 <* v i> + t* 5 < x > ■ - h 
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If the variables are renamed as -follows 


01 “ 


F 


0, = G 


AV : = x 


x = y 


30! 

bav: 


= F 


30 2 

3AVT 


= G. 


x 


^01 

"3T 


= F 


. 3.0, 


2 
3 X 


= G 


Cramers' Rale yields the following solution to C. 4-2: 


6x = 


G F - F G 

y y 

F G - F G 
x y y X 


6(AV£) 


6y = 


F G - G F 
x x 

F G - F G 
x y y x 


= 6(X) 


(C. 4-2) 


Step size control on both components of the step is employed 
as in C. 3. As mentioned in Chapter 4, this simulation is very 
sensitive to the initial values selected and it does not guarantee 
convergence unless the iteration is initiated in the region of the 
stationary point. Good initial values can be selected from the 

•ju 4 

T vs AV' curves obtained from the 1-dimensional iteration C. 3, 
s 

but since these curves supply the information needed, the 
simulation outlined here is only useful as a refinement on data 
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which must he interpolated from the T vs. AV curves. 

C. 5 N- Dimensional Search for Exact Solution 

Given an initial set of values AVp AVg, * * * ^ as 

obtained from the simulation C. 3), an N- dimensional step 
5 (AV^), 6 (AV^), ... 6 (AV.^) can be taken such that the new AV 
array will be driven toward that set of values which satisfy the 
first order necessary condition (3. 2-13) and the velocity constraint 
(3.2-14). In other words, for i = 1, 2, .... (N-l), AV^ must 
satisfy (3. 2-13), and AV^ must satisfy (3. 2-14). Therefore, the 
N constraint equations which must be satisfied are given by 


p i ( AV o- AV 1- AV a) 

= 0 


F 2 Or Av 2- AV s) 

• 

= 0 


• 

♦ 

F N-1 ( AV Nt2‘ AV N-1’ AV n) 

= 0 


G On) 

= 0 

(C. 5-1) 


where AV Q is zero. 

The N- dimensional Newton step 5(AV) that drives these 
constraints to zero is given by the solution to the following set of 
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equations: 


dF, 


c>F 3F 

^ 6(AV 1 )+ Wj 6( ^2 ) + 


, SF, 

• 6(av n } 


= "F, 


SFr, SFr, 

dP 2 = SAV^ 6(A V + 6(AV 2 ) + 


5F ? 

+ 6(av n } 


= -F r 


dF. 


N-l 


^ P N-1 ^N-l ^^TST-1 

^ 6(A V + 3®J- 5(AV 2 )+ --- + ^ 5 < A V = - F N-1 


dG 


tm 6<AV i ) + w; 6 <av 2> + • • • + -^b 6(A V - - G 

1 W (C. 5-2) 


Since G is a function of AV^ only. 


BG 


0 i $ N 


1 i = N ' 


(C. 5-3) 


and since F. is a function of only AV._^, AV^, and AV 


3F 

^£v 7 = 0 For j f (i-1), i, (i+1) 


(C. 5-4) 


The step components 5(A\C) can be solved for from (C. 5-2) 
using Cramer's Rule; but, before giving the result, define the 
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vectors 


x = {av p av 2 , ... av n } 

5x = | : 6(AV 1 ), 6 (aV 2 ), ... 6(aV n )} 

i f 

l = {f v F 2 , ... F N _ r g} (C. 5-5) 

where x is the vector of parameters AV-, 5x is the Newton step 
vector, and if) is the generalized constraint vector. Hence, the 
N- dimensional Newton step is found from the solution to the 
following equation: 


D6x = 


(C. 5-6} 


where D is the matrix of partial derivatives given by 


F 11 

F 

12 

**• F 1N 

F 21 

F 

£ 22 

F 

... r 2N 

F (N-1)1 

F (N-1)2 

• • * F (N-1)N 

0 

0 

... 1 


{C. 5-7) 
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such that 


dip 

[Dl., = 1 


ij 3AV 


The new value for the parameter vector is computed from 


x, ,, .= x, - 5x. 
— k+1 — k — k 


(C. 5-8) 


The Newton step increment is given by Cramer's Rule as 


N . 

.6x, = -^=- 

~ k D 


where each element N. of the n- dimensional numerator vector N 

3 — 

is the determinate of the matrix of partials D with its jth column 
replaced by the constraint vector ij). 


The partial derivatives are computed numerically by 
differencing, and the step size control utilized previously is applied 
to each component of the step increment vector. The analysis 
was terminated when it became evident that a more sophisticated 
step size control was necessary to make the N-dimensional 
constraint converge to within strict tolerances. 
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